TRIPLETS OF GALAXIES: Some Dynamical Aspects 
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, 1. Introduction 

I In celestial mechanics the 3-body problem has a long and rich history, while the 

^ • problem of 3-galaxy systems is rather scarce (see review by Valtonen & Mikkola 

w | 1991). Most of the up-to-date works have addressed the 3-galaxy problem using 

■ a point-like approach, although some explicit-physics simulations have been per- 
I | formed to simulate dynamical friction effects and merging processes (Zheng, et 

al. 1993). These studies have provided important knowledge on the general be- 
haviour of triplets of galaxies that are observed in the sky (Karachentsev 1999), 
and some accordance with observations has been obtained. 
{Nl . However galaxies are not point-like particles, but rather consist of a large 

■ number of stars that in turn can be approximated as point-like particles. Qual- 
itative and quantitative differences result in the dynamics when the 3-galaxy 

' problem is addressed self-consistently; i.e. when galaxies are able to redistribute 

. energy and angular momentum among their stars. Since not using self-consistent 

galaxies casts necessarily some doubts on earlier results, we address here the 3- 
galaxy problem self-consistently and compare some results with observations. A 
^ | full report on this work is in preparation. 

a: 

2. Numerical Experiments 

' Galaxies were modeled after a Plummer sphere with N = 3000 particles each. 

^ , No explicit difference was made as to particles being luminous or dark. The units 

^ ' used here are such that G = M = Rq = 1, where M is the mass of each galaxy 

and i?o its scale-length. To transform iV-body results ('n') to astronomical ones 
('a') we need to choose a real galaxy. We use here a galaxy similar to ours with 
M « 5.5 x 10 11 M Q and R haXo ps 135 kpc (Kuijken & Dubinski 1995, Model B). 
The following transformations follow: 

— ps 13.5 kpc, — ps 5.5 x 10 n M o , — Ps32Myr, — ps 420km/s. (1) 
r n m n t n v n 

The initial positions of the centre-of-mass of galaxies were sampled from a 
homogeneous spherical distribution of radius R max - The evolution of the triplet 
depends, obviously, on the size of this sphere. This radius is taken here as 
an approximate turn-around radius of a density perturbation with the mass 
of a triplet; galaxies are assumed to be already formed. We consider an ini- 
tial virial ratio of 2T/\W\ = 1/4 for this perturbation with velocity dispersion 
a = Vb/2^/3, where V = y/3GM t /5R max and M t = 3M (e.g. Barnes 1985). 
Numerical simulations of galaxy formation tend to indicate that the dark matter 
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Figure 1. Outcomes of a triplet simulation, from initial conditions 
to ~ 20 Gyr in steps of ~ 5 Gyr. (top) A binary and a single galaxy 
results, (bottom) a triple merger is produced. Lines indicate the path 
of the centre of galaxies. Frames are 100 x 100 units. 

background at turn-around had a ~ 20 km/s (Lake & Carlberg 1988). We take 
this as a fiducial value for a in triplets at maximum expansion. Hence for three 
Galaxy-like spirals we obtain i? max ~ 1 Mpc. No common dark matter halo 
has been used in the present simulations. Note that if we increase the mass of 
the triplet by introducing a common dark matter halo, and retain cold initial 
conditions, this will increase i? max proportionally. The collapse time is taken 
here as r co u = n\/ i?j^ ax /2GAf t ~ 794 20 Gyr. These IC's are more appropiate 
for triplets turning around at this epoch. 

We also considered virialized initial conditions (IC's), obtained by assuming 
galaxies to be point-masses. We made 30 simulations and estimated the 1- 
D velocity dispersion a, mean harmonic radius Ru, crossing time t c , and virial 
mass M v (Nolthenius & White 1987), along three orthogonal projections; i.e. 90 
'triplets' are simulated for each of the IC's considered. Simulations lasted ~ 20 
Gyr from turn-around, and energy was conserved for all the runs to ^ 1.5%. 

3. Results 

In Figure |l] we show two typical outcomes of the simulations performed. Qual- 
itatively situations where binaries are formed first, some of them leading to a 
rather quick triple merger (bottom) and others taking a much longer time to 
even form a binary merger (top), are predominantly found. This resembles the 
same instability found in the 3-body problem to form initially a binary. How- 
ever, when self-gravity is considered effects such as a sling-shot are not easy to 
reproduce due to the galaxies' capability to absorb orbital energy. 

In Table [l] we present quantitative results for both types of IC's, taking out 
mergers. Numbers are given in iV-body units. The first row, for each time, are 
average values while in the second median values are given. Times correspond 
respectively to 0.5, 5, 10, 13, 15, and 20 Gyr. As expected, initially virialized 
triplets evolve much slower than cold ones. Collapsing triplets yield a significant 
number of mergers in the time interval (10 — 13) Gyr; i.e. ~ t co \\/2. The average 
or median mass values never overestimate of the true mass of the system; in 
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collapsing triplets. The underestimate in the median is a factor of ^ 3 than 
the true mass during the wide time interval of ~ (5 — 15) Gyr for collapsing 
triplets; the agreement is much better for initially viralized systems. Some 
triplets provided an overestimate in mass along a particular line-of-sight. 

Table 1. Results for Triplets in Isolation 
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Results for the median Rn in collapsing triplets do not agree well with those 
of observations of present day compact triplets (~ 65 kpc). This also happens 
for virialized systems which have a i?H ^ 400 kpc. On other hand, the velocity 
dispersions are always a <^ 50 km/s for both types of IC's, a value which is 
about half of the observed median (~ 100 km/s). Nonetheless, about 10% of 
the simulated triplets have a ^ 100 km/s at t ~ 10 Gyr. Using a larger mass 
and halo size galaxy does not help much in increasing the astronomical median-a 
since velocities scale as v oc y/M/R. 

We recall that when comparing to observational data the assumption that 
Karachentsev's catalog of compact triplets forms an homogeneous and well de- 
fined sample is implicit, but this is not so since e.g. it includes galaxies of different 
luminosities (~ mass) and morphological type. This needs to be considered in 
future studies to make more consistent comparisons with observations. 

On other hand, however, there are some indications that triplets lie in 
the periphery of larger systems of galaxies. In a large-scale structure picture, 
triplets were probably not isolated from tidal perturbations before arriving at 
their present state. Hence it is of interest to estimate probable tidal effects on 
the dynamics of triplets in a Hubble time, even though to a first order. In Fig. |2| 
we present only the results of the cr-distribution under a tidal perturbation from 
a far away 'poor cluster'; the triplet was retained at the same initial distance to 
study the effects of a tidal force. The agreement of the median and average a 
with observations is better when triplets are not considered 'island universes'; 
e.g., (a) ~ 100 at t ~ 10 Gyr when the perturber in Fig. Q is at 3 Mpc. Although 
this value depends obviously on the perturber mass and distance, the results 
manifest the importance external fields can play in the dynamical properties of 
otherwise assumed isolated triplets. They also suggests that the environment 
could have introduced a 'selection effect' in allowing some initially wide triplets 
(i?H ~ 630 kpc) to become compact, and disrupting others, in a Hubble time. 
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Figure 2. Velocity dispersion distribution for collapsing triplets un- 
der the action of tidal perturbation from a point-mass M ~ 10 14 Mq. 
(left) Triplet at 5 Mpc away from the perturber, and (right) at 3 Mpc. 
Dashed-lines are at t = 0, dotted at ~ 10 Gyr, and solid-line at ~ 15 
Gyr. Bars correspond to a « 100,200,300 km/s, respectively. Some 
initially wide triplets get disrupted in t ^ r co ii/2 ~ 10 Gyr. 
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4. Final Comments 

Simulations of isolated triplets show e.g. that in average rather low velocity 
dispersions are obtained when compared to observations. In this scenario an 
underestimate of mass will occur when using the bulk velocity and centre of 
galaxies, probably by a factor of ~ 3. This underestimate can be larger for 
particular triplets if strong signs of interactions are present and if galaxies have 
large dark halos; this situation can be similar for compact groups. Wide triplets 
could have evolved into compact triplets in a Hubble time, although their high 
a remains to be explained. 

Tidal perturbations on the evolution of a triplet, however, might have im- 
portant effects on their dynamics over a Hubble time, and consequently on their 
mass estimation. We suggest that triplet dynamics is closely tied to cosmology, 
and that it is not direct to untangle internal effects from external ones. Obser- 
vational studies that would search e.g. for possible correlations between a and 
the density of the triplets environment could shed light on these issues. 
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